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, For small values of the dimensionless viscosity parameter, namely a ^ 0.1, the dynamics 

' of non-radiating accretion flows is dominated by convection; convection strongly suppresses the 

■ accretion of matter onto the central object and transports a luminosity ^ 10^"^ — 10~^Mc^ from 
, small to large radii in the flow. A fraction of this convective luminosity is likely to be radiated at 

J> ' large radii via thermal bremsstrahlung emission. We show that this leads to a correlation between 

. the frequency of maximal bremsstrahlung emission and the luminosity of the source, i^pcak oc LP'/^ . 

\ Accreting black holes with X-ray luminosities \Q^^LEdd > Lx{0-5 — lOkeV) > 10^'' LEdd are 

■ expected to have hard X-ray spectra, with photon indices F '--^ 2, and sources with Lx < 10^^ Lsdd 
, are expected to have soft spectra, with F 3.5. This is testable with Chandra and XMM. 

o 

2 ' At luminosities less than a few per cent of the Eddington luminosity, black holes can accrete via an 

I advection-dominated accretion flow, or ADAF (Ichimaru 1977; Rees et al. 1982; Narayan & Yi 1994, 1995; 

' Abramowicz et al. 1995; see Kato, Fukue & Mineshige 1998 and Narayan, Mahadevan & Quataert 1999 for 
reviews) . Analytical calculations have shown that such flows should be convectively unstable (Narayan & Yi 
1994, 1995; see also Begelman & Meier 1982); the instability has been confirmed in numerical simulations 
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^ • (Igumenshchev, Ghen & Abramowicz 1996; Igumenshchev & Abramowicz 1999, 2000; Stone, Pringle & 



Begelman 1999; Igumenshchev, Abramowicz & Narayan 2000). 

Igumenshchev & Abramowicz (1999, 2000) found that convection is strong whenever the viscosity pa- 
rameter a is small, roughly a < 0.1. Such "convection-dominated accretion flows," or CDAFs, have a very 
different structure than ADAFs (Stone et al. 1999). The density of the accreting gas varies as p oc R~^/'^ 
(where R is the radius) rather than p oc R^^^^, and the mean radial velocity varies as v oc R~^^^ rather than 
V oc i?^^/^ (Narayan, Igumenshchev & Abramowicz 2000, hereafter NIA). 

NIA showed that, in the numerical simulations, the Reynolds stress due to convection is negative, which 
implies that convection moves angular momentum inwards rather than outwards. This property of convection 
was discussed in the context of thin accretion disks by Ryu & Goodman (1992) and Stone & Balbus (1996). 
In the case of a CDAF, convection is so strong that the angular momentum it transports inwards nearly 
cancels the normal outward transport by viscosity (NIA; Quataert & Gruzinov 2000). As a result, one has 
a nearly static accretion flow in which most of the gas circulates in convective eddies rather than accreting 
onto the central object. For fixed boundary conditions at large radii the mass accretion rate in a CDAF is 
thus much smaller than in a non-convecting ADAF or a Bondi fiow. 
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In this paper we highlight observational consequences of CDAF models, focusing in particular on spectra. 
Some of our results are similar to those of Di Matteo et al. (1999, 2000; DM) and Quataert & Narayan 
(1999a; QN) who followed up the proposal of Blandford & Begelman (1999; see also Narayan & Yi 1994, 
1995) that a significant fraction of the mass in an ADAF would be lost to an outflow/wind, rather than 
accreting onto the central object. The importance of outflows can be parametrized by a radial density profile, 
p oc i?-3/2+P, with p e [0, 1]. The density profile in a CDAF is equivalent to p = 1. Thus DM's and QN's 
spectral models, which considered various values of p, capture many of the features of the CDAF spectra 
calculated here; this is discussed in more detail in §4. 

Convection in CDAFs transports a luminosity Lc ^ 10^^ — lO^^Afc^ from small to large radii; the 
energy is supplied by the small amount of mass accreting onto the black hole. What happens to this energy 
at large radii? NIA and Igumenshchev & Abramowicz (2000) suggested that some, perhaps most, of the 
energy might be radiated from the outer regions of the CDAF as thermal bremsstrahlung emission. We 
present model spectra for such fiows. 

In the next section we present analytical expressions for the luminosity and bremsstrahlung spectrum 
of CDAFs (§2). We then show more detailed numerical calculations (§3.1) and compare CDAF models with 
ADAF models (§3.2). In §4 we summarize and emphasize some implications of our analysis. 



2. Self-Similar Scalings for Bremsstrahlung Luminosity and Spectrum 

In this section we use the self-similar CDAF solution described by Narayan et al. (2000) and Quataert 
& Gruzinov (2000). We employ Schwarzschild units for the radius, i.e. r ~ R/Rs, where Rs — 2GM/c^ = 
2.95 X lO^m cm and m is the black hole mass in solar units. We assume that the accretion fiow extends 
from an outer radius rout down to an inner radius rin = 1 (to model a Schwarzschild black hole) . The sound 
speed of a CDAF is nearly virial, « 0.37c^/r, and the temperature is T = To/r w lO^^K/r; the vertical 
scale height is thus of order the radius, H « 0.6R. 

The density of gas in a CDAF is given by 

/ \ -1/2 

-1/2 I \ 1/2 n\ 

\rout J 

We assume that the CDAF is in steady state with a constant mass accretion rate M. Since p varies as r~^/^, 
the radial velocity must scale as v oc r~^/^. We normalize the velocity such that it is equal to c at 

As noted in §1, convection in CDAFs transports energy from small to large radii. We write the convective 
luminosity as Lc = ecMc^, where Cc ~ lO^'^ — 10"'^ is the convective "efficiency." Let a fraction r/c of this 
energy be radiated at large radii by the gas in the CDAF; this radiation comes out as thermal bremsstrahlung 
emission. The bremsstrahlung luminosity of the CDAF, icDAF, then satisfies 

-^CDAF = VcLc = rjcecMc^. (2) 

For a given Tout; this relation fixes the value of pout B.nd thus M (see below). Equivalently, for a given M, 
equation (|^) uniquely fixes the outer radius of the CDAF, Tout- It also implies an obscrvationally interesting 
correlation between the luminosity and X-ray spectrum of a source. 

The simplest possibility we can consider is that 77c = 1, which corresponds to all the convected energy 
being radiated by the CDAF. This may apply to those objects in which the CDAF is formed by the "evapora- 
tion" of a thin accretion disk, as has been discussed in the context of ADAFs by several authors (e.g., Meyer 
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& Meyer-Hofmeister 1994; Narayan, McClintock, & Yi 1996; Honma 1996; Meyer, Liu, & Meyer-Hofmeister 
2000). If a CDAF is surrounded by a thin disk, and if rjc < 1, the convective luminosity impinging on the 
inner edge of the thin disk would be very much greater than that needed to unbind the disk material. One 
could imagine that would then "eat away" the inner edge of the disk until the condition icDAF ~ is 
roughly satisfied. This condition is not guaranteed, however, since another possibility is that r/c < 1 and the 
excess convective energy is used to eject a substantial fraction of the mass evaporated from the disk. 

When there is no outer disk, it is unlikely that r/c will be close to unity. For example, if a black hole 
accretes spherically from the ISM of a galaxy, as is believed to occur in, e.g., Sgr A* at the center of our 
Galaxy (e.g., Melia 1992; Narayan, Yi, & Mahadevan 1995) or elliptical galaxies at the centers of cooling 
flows (Fabian & Rees 1995; Mahadevan 1997; DM; Quataert & Narayan 1999b), both rout and pout are 
determined by the properties of the host galaxy. It would then require considerable fine-tuning for Lcdaf 
to be equal to L^. It is more likely that we will have ^cdaf < L^., i.e 77^ < 1. The unradiated energy 
{Lc — Lcdaf) will be transported out, heating the external medium and perhaps driving some of the gas 
away in an outflow. 

The bremsstrahlung luminosity of a CDAF is given by 

Ap^T^l^r'^dr, A « 5 x 10^° ergs cm^ g-^ g-i K^^/^ (3) 
Performing the integral and equating icDAF to the convective luminosity, we may solve for p^. 

rjcCcC^ -3/2 10^3 , \ _3 

^° 9-rV2„ ^^T^ilFljgcm . (4) 

2To' RsA mrj^^ ^ ^ 

Substituting this back into the expression for icDAF we find that the bremsstrahlung luminosity of the 
CDAF satisfies 



_ ^CDAF -3/2 f Vcf^c 

/cDAF = -J ~ r^ut [ TWZ^ I > 5 

LEdd ^10 ^/ 

where ^Edd — 1-25 x lO'^^m ergs~^. For rout ^ 1, we see that the luminosity of a CDAF is <C ^Edd- Note 
that equation (^) can also be expressed as a relationship between the accretion rate and outer radius of the 
CDAF: 



^ ^ in -3/2 
~ lOr, 



out {1Q-2) ' 



MEdd VlO 
where MEdd = ^Edd/O-lc^ lO'^^m g is the Eddington accretion rate. 

A straightforward calculation (ignoring the frequency dependence of the Gaunt factor) shows that the 
bremsstrahlung spectrum from a CDAF consists of three parts: 

l/L^ oc hu fcTmin, (7) 

oc 1^-3/2^ kT^in hv kT^^^, (8) 
oc exp(-/ii//fcTmax), ^Tmax < hiy, (9) 

where Tmin ~ T'(''out) is the minimum electron temperature in the flow, obtained at the outer radius, and 
Imax is the maximal electron temperature, obtained close to the black hole. If we include the Gaunt factor, 
the slope is shallower than unity for hv <IC fcTmin. 
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The peak of the bremsstrahlung spectrum occurs at 

/ii^peak ~ ^Tmin = ■ (10) 

^out 

Thus the position of the spectral peak depends on the outer radius. Since the himinosity of the CDAF is also 
related to the outer radius (eq. 5), we have the following relation between the luminosity and the location 
of the peak in the bremsstrahlung spectrum: 

K...=»10(^2f)"°(|S)""^eV. (II) 



3. Detailed Spectra 

The detailed spectra presented in this section were computed using a global dynamical model which includes 
a sonic transition close to the black hole (Quataert, in preparation). The model corresponds to a viscosity 
parameter a = 0.03 and has Cc = 0.0045. The profiles of density and temperature in this model are quite 
close to the self-similar form assumed in the previous section. 

In addition to bremsstrahlung emission, the spectral calculations presented here include synchrotron 
radiation and Compton scattering, assuming thermal electrons. These processes have been included by the 
methods described in Narayan, Barret & McClintock (1997). Cooling from atomic processes has also been 
included, which can be important for T < 10^ K (i.e r > 10^). 

In each model, the electron temperature Tg was computed self-consistently at each radius. Specifically, 

we solved an energy equation as a function of r and ensured that the heating, cooling and energy advection 
of the electrons arc in balance. For the heating, we assumed that a fraction S of the viscous dissipation goes 
directly into the electrons. 



3.1. Models with a Pure CDAF 

Figure 1 shows model spectra corresponding to a CDAF around a W^Mq black hole. All the models 
shown have rjc = 1 and either 5 = 0.01 or S = 0.5. The various curves correspond to different values of the 
outer radius, from rout = 10^ down to 300; the corresponding values of m are given in the caption to the 
figure. (The precise values of m depend on details such as the values of a and Cc, and the global model of the 
gas dynamics, and probably should not be taken too literally. The trends are, however, likely to be robust.) 

The main features discussed in §2 are clearly seen in Fig. 1. Consider first the models with i5 = 0.01, 

shown by dashed lines. We see that the X-ray spectrum is dominated by bremsstrahlung emission and 
consists of two distinct power-law segments on either side of a peak. The location of the peak is highly 
correlated with the luminosity, and follows the scaling given in equation (10). 

At low frequencies, in the radio/submm band, there is a weak secondary peak in the spectrum which 
arises from synchrotron emission by the hot electrons. (This emission was computed assuming that the 
magnetic pressure in the plasma is a tenth that of the gas, cf. Quataert & Narayan 1999c.) The synchrotron 
peak in the models is much weaker than in corresponding ADAF models (Narayan, Mahadevan & Quataert 
1997), where the synchrotron emission is often stronger than the bremsstrahlung emission. This is because 
the synchrotron emission comes from relativistic electrons close to the black hole, whereas the bremsstrahlung 
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emission is from large radii. For a given outer density (and hence bremsstrahlung emission), the electron 
density at small radii in a CDAF is much smaller than that in an ADAF, causing the synchrotron emission 
in the CDAF to be strongly suppressed. This point was made by DM and QN in the context of winds from 
ADAFs, which also cause a reduction in the central density. 

For higher values of M (or luminosity), Compton-scattering of synchrotron photons becomes important; 
this fills in the valley between the synchrotron and bremsstrahlung peaks. 

The ions in the above models are close to virial: Tj ~ 10^^ K near the black hole. The electrons are 

comparatively cool, with a maximum temperature T^^max < 10^° K. Models with S = 0.5 (solid lines in Fig. 
1) have much stronger electron heating and therefore significantly hotter electrons: ^e,max ^ 

10" K. This 

causes the synchrotron peak to become much more pronounced compared to the 6 = 0.01 models, though 
the emission is still less than in an equivalent ADAF model. The increased synchrotron emission and the 
hotter electrons both lead to an enhanced contribution from Compton scattering. Nevertheless, the peak of 
the spectrum is still dominated by bremsstrahlung emission, and so the scaling of the position of the peak 
with luminosity survives. 

Figure 2 shows spectra for 6 = 0.01 and 6 = 0.5 for several values of the convective efficiency, r]c- the two 
sets of models correspond to m = 10~^ (solid hues) and 10~^ (dashed lines). In all the model sequences the 
synchrotron emission in the radio is relatively unaffected by changes in rjc] this is because the synchrotron 
radiation arises primarily at small radii and depends only on rh and 6. For the large m model, varying rjc also 
has little effect on the X-ray spectrum. This is because most of the X-ray emission arises from small radii 
via inverse Compton scattering, rather than from large radii via bremsstrahlung. By contrast, for smaller rh 
(and/or smaller 6), inverse Compton emission is less important. The X-ray emission is then dominated by 
bremsstrahlung, which decreases with decreasing rjc- 

Figure 3 shows the predicted X-ray photon index F as a function of X-ray luminosity in Eddington 
units for 6 = 0.5. Four sets of models are shown: r/c = 1, 0.5, 0.2, 0.1. The 2-5 keV spectral indices were 
calculated by comparing the model luminosities at the two ends of the range, 2 keV and 5 keV, respectively. 
Note the large increase of F with decreasing luminosity. This arises because the peak in the spectrum moves 
to lower energies with decreasing luminosity (cf. eq. 10). Because the value of rjc is unknown, there is 
some uncertainty regarding the value of the luminosity at which the transition in the value of F occurs . 
For the range of rjc considered here, covering an order of magnitude, we find that for Lx/Lsdd 10~^ the 
X-ray spectrum is hard, with a photon index F ~ 2, while for luminosities Lx/LEdd ^ 10~® the spectrum is 
significantly softer, with F 3.5. 

3.2. CDAF vs. ADAF 

Figure 4 shows a comparison of CDAF and ADAF spectral models for S = 0.01 and d = 0.5. The 
accretion rates in the models are adjusted so that the 1 keV X-ray luminosities are equal to either 10^'' or 
10*° erg s~^. In order to uniquely fix the models, we set rjc = 1 m the CDAF models and choose the outer 
radii of the ADAF modles so that local viscous dissipation balances radiative cooling at this radius [f <^l/2 
in the notation of Narayan & Yi 1994). 

For small 5, Figure 4 shows that the CDAF and ADAF spectra are quite different. The X-ray emission 
in the CDAF is dominated by bremsstrahlung, while the emission in the ADAF is dominated by inverse 
Compton emission. In addition, the CDAF has a much lower synchrotron flux in the radio, because the 
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density of gas close to the black hole is much smaller. 

For very low luminosities (< 10~^LEdd), the difference between the CDAF and ADAF models at large 
S is similar to that at small S. By contrast, for luminosities > 10~^LEdd the CDAF and ADAF spectra arc 
very similar if S is large. This is because the entire spectrum is produced over a small range of radii close to 
the black hole. Fixing the 1 keV X-ray luminosity determines the properties of the flow in this region and 
thus the entire spectrum is rather similar. 

4. Discussion 

The main results of this paper are given in equations (6)-(10) and in the Figures. Equation (10) shows 
that a CDAF has a peak in vLi, in the X-ray band, whose position is correlated with the Eddington-scaled 
luminosity. Correspondingly, the spectral index of the X-ray spectrum varies with luminosity. At high 
luminosities, the peak in the spectrum is at higher energies than the typical band in which observations 
are carried out (few keV), and the spectrum is hard with a photon index F of order 2 (Fig. 3). At lower 
luminosities, the peak shifts to lower energies, and the observed band would correspond to energies above 
the peak. The spectrum then becomes very soft, with F taking values of order 3.5. 

There is some uncertainty in the quantitative details predicted here because the values of two parameters 
are not well determined. The parameter ec measures the fraction of the rest mass energy of the accreting 
gas that is carried outward by convection: = CcMc^. Our calculations correspond to eg = 0.0045, as 
determined from a global CDAF model with viscosity parameter a = 0.03 (Quataert, in preparation). Other 
values of a would give other estimates of Cc, though the uncertainty in Cc is probably no more than a factor 
of a few either way. 

A more serious uncertainty is in the parameter i]c, which measures the fraction of the convected lu- 
minosity Lc that is radiated by the CDAF at large radii as thermal bremsstrahlung emission: Lcdaf = 
r/c-t'c = r]c^cM(?. Figure 2 shows results spanning a range of r/c from 0.01 to 1. If r]c is as small as 0.01 or 
even smaller, an accreting black hole would have a very small X-ray luminosity from thermal bremsstrahlung 
emission and would be extremely difficult to observe. This might explain highly underluminous galactic nu- 
clei such as that in our own Galaxy. For models with low ric, a relevant question is what happens to the part 
of the energy convected outwards that is not radiated, ([1 — r?c]ic ~ L^)- If this energy is transported out to 
an external medium and radiated there, that contribution may dominate the observed spectrum (assuming 
the beam of the telescope is larger than the size of the emitting region). The luminosity of the source would 
then be larger than we have predicted and the spectrum would be softer. 

The spectral models of CDAFs presented here are qualitatively similar to the ADIOS (ADAF -|- wind) 
models presented by DM and QN. The relative prominence of bremsstrahlung emission in the X-ray band and 
the suppression of synchrotron emission in the radio (as suggested by observations; see DM) are characteristic 
of density profiles flatter than the ADAF scahng of p oc r~^/^ (be they CDAF or ADIOS). In fact, the 
qualitative similarity of ADIOS and CDAF spectra implies that direct signatures of outflowing gas (and 
measrements of outflow mass rates) would be needed to confirm the ADIOS model. The primary difference 
between the CDAF and ADIOS models is that we have a unique radial density profile [p oc r~^/'^) in a 
CDAF, instead of a family of models {p oc r~^/'^~^P with < p < 1) in an ADIOS; also, the energy which 
powers bremsstrahlung emission at large radii is provided by convective transport from small radii in a 
CDAF, rather than by local viscous dissipation in an ADIOS. 
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All the models presented in this paper correspond to a two-temperature plasma. For the models with 
S = 0.5, however, the temperatures of the ions and electrons are similar, differing by less than a factor of 
10. We have computed models with a one-temperature plasma, assuming that some plasma process other 
than Coulomb collisions rapidly equilibrates the ion and electron temperatures (cf Begelman & Chiueh 
1988). The results arc not very different from those obtained with two-tcmpcraturc models with S ~ 0.5. 
Comptonization does, however, become more important relative to bremsstrahlung, especially for low values 
of r]c- In addition, the synchrotron emission in the radio is somewhat more prominent. 

An important assumption of our analysis is that the electrons are thermal. Because of the low density in 
a CDAF, thermalization through Coulomb collisions and synchrotron self-absorption is virtually non-existent 
(cf. Mahadcvan & Quatacrt 1997); it is therefore possible to retain a power- law distribution of electrons. 
Moreover, such nonthermal acceleration is expected in the coUisionless magnetized plasmas present in CDAFs. 
A power law tail of electrons would significantly modify the predicted synchrotron spectrum (e.g., Mahadevan 
1998; Ozel, Psaltis, & Narayan 2000). In addition, the power-law electrons could contribute to X-ray emission 
through cither synchrotron or nonthermal inverse Compton emission. Depending on the imccrtain efficiency 
of electron acceleration, this nonthermal emission could be more important than bremsstrahlung emission, 
particularly for small values of r/c- 

Models in which thin accretion disks "evaporate" to form hot inner flows have received considerable 
attention in the context of dwarf novae (Meyer & Meyer- Hofmeister 1994), soft X-ray transients in quiescence 

(Narayan, McClintock, & Yi 1996; Honma 1996; Meyer, Liu, & Meycr-Hofmeistcr 2000, Rozanska & Czerny 
2000), and low luminosity AGN (Lasota et al. 1996; Quataert et al. 1999). The energy for "evaporation" is 
ultimately thought to originate in the hot inner flow. Convection in CDAFs transports a significant amount 
of energy to large radii. CDAFs should thus be much more efficient than ADAFs at evaporating the outer 
disk. The implications of convective energy transport for theoretical models of the transition from thin disks 
to hot flows deserve further investigation. 
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Fig. 1. — Model CDAF spectra for a W^Mq black hole, taking r]c = 1. Dashed lines are for a fraction S = 0.01 

of the viscous energy heating the electrons. From top to bottom, the curves correspond to models with 
(log rout, log m) = (2.5,-3.29), (3.0,-3.86), (3.5,-4.49), (4.0,-5.19), (4.5,-5.99), (5.0,-7.15), (5.5,-8.36), 
(6.0,-9.70), respectively. Solid hues are for 6 = 0.5. From top to bottom, the curves correspond to 
models with (log rout , log m) = (2.5,-3.38), (3.0,-3.91), (3.5,-4.50), (4.0,-5.19), (4.5,-5.99), (5.0,-7.15), 
(5.5, -8.36), (6.0, -9.70), respectively. 



- 11 - 




Fig. 2.— Model CDAF spectra for a lO^M© black hole. Left panel: Models with 5 = 0.01. The solid 
curves correspond to m = 10"'* and, from top to bottom, (logrout, ?7c) = (3.1,0.96), (2.9,0.55), (2.7,0.31), 
(2.5,0.16), (2.3,0.074), (2.1,0.019), respectively. Dashed curves correspond to to = 10"® and, from top to 
bottom, (log rout, r?c) = (4.5,0.97), (4.3,0.45), (4.1,0.22), (3.9,0.11), (3.7,0.056), (3.5,0.029), (3.3,0.015), 
(3.1,0.0080), respectively. Right panel: Models with 5 — 0.5. The solid curves correspond to m = 10~^ and, 
from top to bottom, (log Tout, ??c) = (3.0,0.86), (2.8,0.54), (2.6,0.34), (2.4,0.20), (2.2,0.094), respectively. 
Dashed curves correspond to to = 10"® and, from top to bottom, (logrout, ??c) = (4.5,0.97), (4.3,0.45), 
(4.1,0.22), (3.9,0.11), (3.7,0.059), (3.5,0.032), (3.3,0.018), (3.1,0.011), respectively. 
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Fig. 3. X-ray photon index vs. X-ray luminosity for models with 6 = 0.5. Solid, dotted, dashed, and 
long-dashed curves correspond to rjc = 1, 0.5, 0.2, and 0.1, respectively. 
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Fig. 4. — Comparison of ADAF and CDAF model spectra; rjc = 1 for the CDAF models. The accretion rates 
in the models have been adjusted so that the 1 keV luminosities are equal to either 10'^'^ erg s"-'^ or 10^" erg s^^ 
(solid circles). Left panel: Models with 6 = 0.01. The ADAF models (solid curves) correspond, from above, 
to (log Tout) log rh) = (4.9,-2.74) and (5.2,-3.64), respectively, and the CDAF models (dashed curves) cor- 
respond to (log Tout , log m) = (2.9,3.75) and (4.7,-6.37), respectively. Right panel: Models with 6 = 0.5. 
The ADAF models (solid curves) correspond to (log rout , log m) — (5.5,-3.53) and (5.6,-4.20), respec- 
tively, and the CDAF models (dashed curves) correspond to (log rout) log to) = (3.3, —4.26) and (4.7, —6.37), 
respectively. 



